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Abstract 

We solve exactly the equations of motion for linearized gravity in the Randall- 
Sundrum model with matter on the branes and calculate the Newtonian limit in it. 
The result contains contributions of the radion and of the massive modes, which change 
considerably Newton's law at small distances. The effects of "shadow" matter, which 
lives on the other brane, are considered and compared with those of ordinary matter 
for both positive and negative tension branes. We also calculate light deflection and 
Newton's law in the zero mode approximation and explicitly distinguish the contribu- 
tion of the radion field. 

Keywords: Kaluza-Klein theories; branes; linearized gravity; Newtonian limit 

1 Introduction 

The present-day approach to Kaluza-Klein theories is based on the idea due to Rubakov and 
Shaposhnikov of localization of fields on a domain wall |]l], ^ . In a first step it is natural to 
drop the mechanism of localization of fields and to treat the domain wall as an infinitely thin 
object, i.e. as a membrane, and to consider the effects due to the gravitational interaction 
of such 3-branes. 

A particular realization of this scenario was proposed in paper [Q] . In this paper an exact 
solution for a system of two branes interacting with gravity in a five- dimensional space-time 
was found. This model is called the Randall-Sundrum model (usually abbreviated as RSI 
model), and it is widely discussed in the hterature (see Refs. 0, || for reviews and references). 
A consistent analysis of this model without matter on the branes was made in . However, 
the equations of motion for metric fluctuations, when there is matter located on the branes, 
have not been studied in detail yet. In paper [0] the Randall-Sundrum model with matter 
on the branes was discussed for the cases of both one and two branes (the former is usually 
called RS2 model). This paper utilized the Gaussian normal coordinates and the bent-brane 
formalism, which was argued to be inconsistent in RS2 model in paper p|. Moreover, the 
use of Gaussian normal coordinates mixed the contributions of the graviton and the radion 



fields to the four-dimensional gravitational field. In the present paper we will decouple the 
equations for the graviton and the radion fields and solve exactly the equations of motion 
for the RSI model with matter on the branes. Then we will calculate the Newtonian limit 
in this model and the light deflection by a point-like mass in the zero mode approximation. 

The Randall- Sundrum model describes the gravity in a five- dimensional space-time 
E with two branes embedded into it. We denote the coordinates by {x*^} = {x'^,?/}, 
M = 0, 1, 2, 3, 4, /i = 0, 1, 2, 3, the coordinate x"^ = y parameterizing the fifth dimension. It 
forms the orbifold S^/Z2, which is realised as the circle of the circumference 2R with points 
y and —y identified. Correspondingly, we have the usual periodicity condition in space-time 
E, which identifies points {x,y) and {x,y + 2nR), and the metric qmn satisfies the orbifold 
symmetry conditions 

9^u{x,-y) = gf,u{x,y), (1) 
gi^iix, -y) = -gf^iix,y), 
5-44(2;, -y) = gAi{,x,y). 

The branes are located at the fixed points of the orbifold, y = and y = R. 
The action of the model is 

S = Sg + S, + S2, (2) 

where Sg, Si and ^2 are given by 

Sg = / (i? - A) d'xdy, (3) 

IottG Je 



51 = Vi y/-g6{y)d xdy, 

Je 

52 = V2 y^S{y - R)d'^xdy. 



Here (jf^u is the induced metric on the branes and the subscripts 1 and 2 label the branes. 
We also note that the signature of the metric gMN is chosen to be (—1, 1, 1, 1, 1). 
The Randall- Sundrum solution for the metric is given by 

ds^ = gMNdx^dx^ = e^'^'^y^r^^M^dx'' + dy^, (4) 

where rj^^ is the Minkowski metric and the function (T{y) = —k\y\ in the interval —R < y < R. 
The parameter k is positive and has the dimension of mass, the parameters A and Vi 2 are 
related to it as follows: 

A = -12k\ Vi = -V2 = 

AttG 

We see that brane 1 has a positive energy density, whereas brane 2 has a negative one. The 
function a has the properties 

d,a = -k signiy), ^ = -2A;(5(t/) - 5(y - R)) = -2kS. (5) 
Here and in the sequel 84 = 
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We denote k = where G is the five- dimensional gravitational constant, and 

parameterize the metric qmn as 

Qmn = iMN + khuN, (6) 

Hmn being the metric fiuctuations. Substituting this parameterization into (0) and retaining 
the terms of the zeroth order in k, we get the second variation action of this model 0. It is 
invariant under the gauge transformations 

^M7v(a;, y) = hMNix, y) - (Vm^7v(x, y) + VN^uix, y)), (7) 

where Vm is the covariant derivative with respect to the background metric •Jmn, and the 
functions ^N{x,y) satisfy the orbifold symmetry conditions 

e{x,-y) = eix,y), (8) 
^^{x,-y) = -^^{x,y). 

With the help of these gauge transformations we can impose the gauge 

h^4: = 0, /l44 = h44{x) = (9) 

which will be called the unitary gauge (see |^). We would like to emphasize once again that 
the branes remain straight in this gauge, i.e. we do not use the bent-brane formulation, 
which allegedly destroys the structure of the model (this problem was discussed in P). 
The general form of interaction with matter is standard. 



^2dx, (10) 



where T^^ and T^j^ are energy-momentum tensors of the matter on brane 1 and brane 2 
respectively: 

rpl,2 _ 2 _ l,2rl,2 

1^" ~ S^f^i^ ^l^" ' 



As follows from formula ([T0|) , /i^j, is the only physical field of the model, since only this 
field interacts with matter on the branes. Obviously, the unitary gauge conditions (^ do 
not fix the gauge of this field. In fact, after imposing these gauge conditions there remain 
gauge transformations of the form 

= e'^e.ix), (11) 

which change the longitudinal components of the field h^^. Nevertheless, it turns out that 
it is convenient to solve the equations of motion for linearized gravity in the unitary gauge 
and then to choose an appropriate gauge in our four- dimensional world on the brane. We 
will use the de Bonder gauge for the field h^^i, on the branes, which corresponds to the choice 
of harmonic coordinates. 
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2 Newtonian limit 



The equations of motion for different components of the metric fluctuations in the unitary 
gauge take the form (see 0): 
1) /iz/-component 

\ [Opd^h^, - d.dPh^, - d^d^hp, + - (12) 



2) /i4- component, 

d^id^h - 9" V) - 394^5^0 = 0, (13) 

which plays the role of a constraint, 

3) 44- component 

^{d'^d'^h^, - d^d^h) - ^d^adji + 6A;V = 0, (14) 

with T^y being the energy-momentum tensor of the matter and h = '^^^h^^,. In what follows, 
we will also use an auxiliary equation, which is obtained by multiplying the equation for 
44-component by 2 and subtracting it from the contracted equation for /iz/-component. This 
equation contains h and only and has the form: 

2 + 2diadih - 8A;V + 8k^S + d^d^^cj) = '^T'^. (15) 

3 

If Tfj^i, = 0, the physical degrees of freedom of the model can be extracted by the substi- 
tution IBI 



h^u = h^u + 7^^(cr ~ c)0 + ^ (^cr - c + - + -e ^'"J d^d^cf). (16) 

with c = ^ikn_i - It turns out that the field b^uix^^y) describes the massless graviton |^ 
and massive Kaluza-Klein spin-2 fields, whereas 0(x) describes a scalar field called the radion. 
Apparentely, the radion field was first identified in Ref. ]TD|] (see also |TT|) and discussed in 

El 0,0,0- 

However, the situation is rather different, when there is matter on the branes. Let us 
first consider the case, where matter is located on the brane at 0, i.e. the energy-momentum 
tensor is of the form T^,^ = t^u{x)5{ii). The substitution, which allows one to decouple the 
equations, looks like 
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After this substitution equations (|T3|), (JT^) and (|T^), rewritten in the flat metric (i.e. u 
rj'^'^u^^), take the form: 



d,ie-'-id''u^,-d^u)) = 0, 



(19) 



d,ie"'d,ie-''^u)) + 2^5{y) [e^'^^ - l] 00 = ^t^^ 5{y). (20) 



Let us consider Fourier expansion of all terms of equation ( pO]) with respect to coordinate 
y. Since the term with the derivative ^4 has no zero mode, this equation implies that 

□0 = ^t, t^t;:, (21) 

di{e-'^''u) = 0. (22) 
From the last equation and equation ([T^) it follows that 



d'^u^, = e2M^(x), (23) 
d^u: = e'"'B^{x), (24) 

where A^(x) and B^{x) depend on four- dimensional coordinates only. It is easy to see that 
the remaining gauge transformations (^) allow us to impose the de Bonder gauge condition 
on the field 

(u^u - ^V^uu^ = 0. (25) 

Having imposed this gauge, we are still left with residual gauge transformation 

= e'"e^(a;), De^ = 0. (26) 

The gauge transformations with satisfying these conditions are important for determining 
the number of degrees of freedom of the massless graviton. It follows from equation (p!9D 
that 

Now let us consider /iz/-equation. After substitution ([T7|) and in the de Bonder gauge it 
takes the form 

^ e~'^''Uu^^ + l:d4,d4Uf,^ - 2k'^u^^ - did^au^^ = (28) 



2 ^ 2 



^jfiu^yyj . 12(1 - e-2^-f^) V □ 

where 

r,,„ = r,,„ I n,„, ^ 



f,. = t,.-l(v,u-^]t. (29) 



One can see that is transverse-traceless. Here and below the inverse d'Alembertian is an 
integral operator uniquely defined by the radiation conditions. 
To solve this equation, we make the following substitution 



Vnu + 



kk 



which does not violate (BSl) that takes the form 



d^'v^, = 0, 



□ 



(30) 



(31) 



Substituting (|^) into (pSj) , we get 



1 92 



2dy 



2 "^^J■>^ 



dy 



--fi^Jiy) 



(32) 



This equation can be solved exactly. To this end, let us first perform Fourier transform of 
the field v^u with respect to the a;-coordinates 



1 -2., 1 



2dy^ 



where = —pi + . 



dy 



-'^UuS{y), 



First let us solve this equation in the bulk (see 
3f) looks like 



ok\y\ 



k 



(33) 

Here the solution of equation 

(35) 



k 



,p'>0 



Substituting it into equation (|33D and comparing the terms at the boundaries, we get the 
values of the constant tensors and -D^jy. Having got them, we obtain the solutions of the 
/xi/-equation with matter on the positive tension brane: 
l)p2 < 



X 



'V,^u{p,y) 



p2 



J2 



-JLek\y\ N, 



(36) 



Ji 



Ji 



'P^ f.kR 
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2)p2 > 



1 / _ Pfip,. 



K 



(37) 



Since we want to calculate the Newtonian limit, t^j^ is proportional to 6{po) in this case. It 
means, that we need the solution for p"^ > 0. 

When matter is on brane 2 {a.t y = R), all the reasonings are the same, as presented 
above. The full substitution looks like 



1 



1 



(38) 



kk 



dadu 



□ 



t. 



The gauge conditions are the same too. The equations of motion for the field and solution 
for the field v^i, with matter on brane 2 can be derived analogously and read 



□c 



K 



:t. 



l)p2 < 



Vf^u — I t 



k 1 



J 



,2 pi'R 



X 



J2 



Ji 



Ji 



,„2 



2)p2 > 



(39) 



(40) 



J 2Vp^e 



,2 p'^^ 



X 



A(#)A-,(#e-)-/,(#e«)A-,(# 



(41) 



An important point is that these equations are written in the coordinates {x^}, which are 
Galilean on brane 1 (not on brane 2) and are inappropriate for studying physical effects on 
brane 2 (we recall that coordinates are called Galilean, if g^^ = diag{—l, 1, 1, 1) [|16| ). 
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Obviously, when there is matter on both branes, the solution for the metric fluctuation 
is just the sum of solutions for each brane separately, which follows from the linearity of the 
equations. 

Now let us examine the four- dimensional theory on the branes. First we consider the 
case when matter and observer are located on brane 1. It is easy to see that h^^ (p^Tl), (|30D 
does not satisfy the de Bonder gauge condition. The residual gauge transformations (^6|) 
are not sufficient to pass to this gauge. But since we consider only the effective theory on 
brane 1, we can drop the De^ = condition, which fixes the gauge for the field Vfj_iy. Then 
we can pass to the de Bonder gauge condition for the field h^^, on brane 1 with the help of 
the gauge functions of the form 

1 /^^ 1 \ a , kR da ^ 



Having been made, these transformations result in 



dnd. 



hi^Ay=o = lj/=o + VtiukR(j) + 2kR—^(f). (43) 

This formula gives the four- dimensional gravitational field on brane 1, and it takes into 
account the contributions of the massless and of the massive gravitons and the contribution 
of the radion. 

Now let us calculate the Newtonian limit in this case. Let us consider a static point mass 
with energy-momentum tensor too = M6{x), tofc = Uj = (and too = 27rM5(j9o)), M denoting 
the inertial mass. We need to calculate only /loo-component of the metric fiuctuations (see 



and 



f,„(p, 0) = = — ^ — ^ ) ( i ^ , ^ , X (44) 

Air 

x—M6{po) 



1 k M , , 

hoo{x,0) = voo{x,0) -kR— — — , 45 

Att bR r 

(see (|4^) ) where r = \^ and 

^oo(x, 0) = ^ J e-'^--P'^^Vooip, 0)d^p. (46) 

Using the last relation and equation ( ^4]) one can easy find that 

1 47r An 



t;oo(x,0) = 777^ — —M / g{p)sin{pr)pdp, (47) 



9{p) 



2p/i(f)^i(fe'=^)-/i(fe'=«)i^i(f) 
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where p = yfP- It is impossible to evaluate the integral in (|47| ) analytically. But we can 
estimate the integrand in the following cases: 1) p « e~^^k, 2) e'^'^k << p << k and 
3) k « p. Using these estimates in the intervals < p < e~^^k, e~^^k < p < k and 
k < p < oo respectively, we can estimate the integral in equation (^). 

First, utilizing the recurrent relations for the McDonald functions, we represent the factor 
g{p) in the integrand in the following form 

. . ^kk_ k Ko (f ) /, (f e^«) + Jo (I) K, (f e^^) 
'^'^ V' 2p/,(|)ir,(|e-)-/,(|e-)K,(f)- 

With the help of this transformation we have picked out the contribution of the zero mode 
to vqq. Substituting asymptotic formulas for McDonald functions into (^81) , we get 



9{P) ~ ^ + 



lik 1 k In (P\ -n^^p-^^h 

^2kR_ip2 2k[l-e-2kRy"' \k) ' P ^ 1^ 

-f/n(f) , e-^^k«p«k (49) 

, k «p. 



2p 



All the integrals arising after substituting (^9]) into (|47|) can be calculated analytically. The 
integrals of the form sin{pr)dp are calculated with the use of the regularization 

f°° f f°° _ \ cos(ar) 
j sin{pr) dp = lim I J e ^^sm{pr)dp \ = . (50) 

It is well known that the gravitational potential is expressed through the ^foo- component of 
the metric as goo = — 1 — 2V ll^. Since goo = — 1 + khoo, we get 



2 

Thus, using equations (^31), ( ^71 ) and (^), we get 



V = -'^hoo. (51) 



y - -T^^ (l -H ^ - 3^ I'Mtr) - s,ikr)]) , (52) 



where Gi = and 

sint 



i-e- 

si{h) = I dt. 

Jo t 

Thus, we have examined the case of the mass and the observer being located on brane 1. 
But there are three more possible cases to be examined. It is the case of " shadow" matter, 
when the mass is located on brane 2 and the observer is located on brane 1, the case of the 
mass and the observer being located on brane 2, and the case of the mass being located on 
brane 1, whereas the observer being located on brane 2 ("shadow" matter effect as well). 
All the calculations in these cases are the same, as described above. But if the observer 
is located on brane 2, it is necessary to pass to Galilean coordinates on brane 2 to get a 
correct result. This problem was discussed in detail in paper |p. The energy-momentum 
tensor t^j, takes the form too = M6{x), tok = Uj = in the Galilean coordinates on the brane 
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the mass is located on. In the formulas presented below (in the Galilean coordinates on the 
observer's brane) the energy-momentum tensor always has this canonical form, and passing 
to the Galilean coordinates is taken into account by a factor in front of the brackets. 

Now let us discuss, how to pass to Galilean coordinates on brane 2 by the example of the 
case when the mass is located on brane 1 and the observer is located on brane 2 ("shadow" 
matter effect). We consider the following form of the Fourier transform with the integrand 
appropriate for this case 



V(a^o,x) = J^,B,, I e'^-^^'^' g{p)6{po)d'p, (53) 



with B^i, being a constant. One can easily find, that in Galilean coordinates on brane 2 this 
expression looks like 

A'^^izo, z) = j^,B,^ e-''^ I e^^^'^^''^'g{e-'\)5{qo)d\ (54) 

where {z'^} are Galilean coordinates on brane 2. In particular, this means, that the intervals, 
in which we have to estimate the integrand, change. Now we have to choose the following 
domains of estimate: q << k, k << q << e^^k and e^^k « q. 

Another difficulty, which arises in this case, is that one of the integrals cannot be evaluated 
analytically. It has the following form 

/ e~fc A/gsin(gr) = A;^/^ / y/ism{tkr)e~^ dt, 
Jk Ji 

where r = But we can estimate it in the limiting cases 1 << kr and kr << 1. 
1) 1 << kr. In this case (|55|) can be estimated by integration by parts 

_^ cos(/cr) 



(55) 



/oo 
Vt sm{tkr)e~^ dt e 



(56) 



kr 

2) fcr << 1. In this case we can make the expansion of sin{tkr) for the small argument 



/oo POO 
Visin{tkr)e^^ dt ^ / \/isin{tkr)e^* dt - 



(57) 



Vt (^krt - i^LLL + ... ] dt. 



The first integral was taken from the table of integrals |jT8| and the second one was calculated 
numerically. The result reads as follows 



(5i 



And the last very interesting point of "shadow" matter cases is passing to the de Bonder 
gauge. Let us show it in the case of the field h^y\y^^. It does not satisfy the de Bonder 
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gauge condition on brane 2 (see (|38D ), analogously to the case of matter located on brane 1. 
With the help of the gauge functions of the following form 

^2kR Kp 2kR fi 

we can pass to the de Bonder gauge. But after this transformation h^i^l^^j^ appears to have 
the following form 

Vlj;=fi = ^M^'lj;=i?- (60) 

This means, that h^u\y^j^ satisfies transverse-traceless gauge conditon 

d''h,.\^^^ = 0, (61) 

K\y=R = 0, 

which is quite natural, because there is no matter on brane 2. 
Thus, the Newtonian limit in this case has the following form 

1) 1 << kr 

2) kr « 1 

V . _G.e-M. A (i , ^ _ 0,^ , 0.099^i(..)=) . (63) 

^ -Gie-^^Mk (2.05 + Qm{krf) . 

Below the results for the last two cases are presented. All calculations were made similarly 
to the presented above. 

I) The mass and the observer are located on brane 2: 



V ^ -G2— 1 - e'""" + — — si{kr 



r \ Sir Svr 



cos(A;r) sin(/cr) 



+ 



kr {kry 



(64) 



where G2 = ^2kR_i is the effective gravitational constant on brane 2. 

II) The mass is located on brane 2 and the observer is located on brane 1 ("shadow" 
matter also): 

1) e''^ « kr 

V « -le-G,^ f 1 + ,/|e-e-H5?(^^ | . (66) 



2) kr « e 



kR 



V . (l + ^ - 0,2^1+ 0.099e—y|i(*r)" ) « (66) 

^ -G2Mk (2.05 + 0.02 e-^^\krf) . 
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We would like to note that in all formulas for the Newtonian limit, presented above, param- 
eter M denotes the physical mass in Galilean coordinates on the brane the mass is located 
on. 

In paper |]^ possible values of parameter k were discussed. It was shown that the value 
k ~ ITev is the only one admissible for a correct physical interpretation of the theory on 
brane 2. This means, that for r ~ 1cm the value kr ~ 10^^ and the terms proportional 
to sin(A;r) and cos{kr) can be dropped as negligible. For example, since si{t — > oo) — > |, 
equation (|4|) takes the form 

V ^ -G,y (l + le^^^) . (67) 

The same arguments can be applied to other cases of matter and observer disposition. Since 
the contribution of massive modes are negligible for r » 10~^ ''cm, only the zero modes 
constitute the Newtonian limit. Unfortunately, in formulas (p2D, (|62|), (|6^) and (^) the 
contributions of the massless graviton and of the radion are " mixed" . This problem can be 
solved with the help of the zero mode approximation. 



3 Zero mode approximation 

Let us consider equations (|12]), ([T3| ) and ([13|) with the energy-momentum tensor of the form 



T^u = t^y{x)6{y) (in the unitary gauge). As described above, with the help of (p!7D these 
equation (except for the /iz/-equation) take the form (p!8D , (|T9| ) and (^Ol). Then we can get 



equations ( |2T| ) and ( p2D as it was made above. It is well known that the field u^j^y in the 
presence of matter is a combination of zero and massive modes 0, whose eigenfunctions are 
orthogonal ||^. In particular, the zero mode can be represented as = e^^a^y, where a^i, 
depends on x only. It also means that with the help of the residual gauge transformations 
= e^"e^{x) it is possible to impose the gauge condition 



d'^u-;, = 0, (6^ 

on the massive modes u^^^ and the de Bonder gauge condition on the zero mode 



( - i?7^^a ) = 0. (69) 



After imposing this gauge, we are still left with residual gauge transformation (|26|). It follows 
from equations (|19D, (|l]),(|8D, (|2|) that 

Ua = — ^t. (70) 

Now let us consider equation (p!2D . Making substitution (|T7D with condition (plj) and 
passing to the gauge (pSF), (pUD we get 



^□(«,. - lv,ua) + ^e-2-n«- + ^d,d,u^^ - 2k'u^^ - d,d,au^^ = (71) 
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Since we are going to calculate the Newtonian limit and light deflection in the zero mode 
approximation, we have to find an equation for the field a^^. If we multiply equation ( |7TD 
by e^'^, integrate it over y and take into account the orthonormality condition for the wave 
functions of the modes, we get 

0{a^u - ^ri^ua) = - (^_g-2fcR) V- (72) 
The equation for massive modes takes the form 



le-2-n«- + \d,d,u^^ - 2k\^^ - d,d,au^, = -^t,J{y) + (73) 
kke^^^ kke^^^ ( d.A, \ k ( d„d. 



^2(e2^«-l)^^"^ 6(e2fc«-l) (v ^ "If^-J^^iy)- 

We note that equation (|7TD has been solved exactly in Section 2. But it turned out to 
be impossible to evaluate analytically all the arising integrals even for a simple form of t^i/ 
(for example, for a static point mass). 

When matter is on brane 2 {at y = R), all the reasonings are the same, as presented 
above. The substitution has the form (^). The gauge conditions are the same as in case 
of matter on brane 1 too. The equations of motion for the fields a^,^, in the presence of 
matter on brane 2 can be derived analogously and read 

loia,,-\v,ua) + \e~^''^u^u + \d,d,u^,-2k^u';:,-d,d,au^,= (74) 

r/ 7-.\ ( '^U.'^y \ CI 

= --t^J{y-R) + g(g2fcR _ 1) - ^M'- )K- 

- ^(^-,,.)t;:%-/2), 

1 hhk 

□0 = —f"^. (76) 
^ 6R ^ V ; 

We would like to note once again that these equations are written in the coordinates {x'^}, 
which are Galilean on brane 1 and are inappropriate for studying physical effects on brane 2. 

Now we are ready to find the Newtonian limit and light defiection by a point-like static 
mass. Let us first make it for brane 1. The substitution (0) in the zero mode approximation 
is 

= - g2fcfl^ I ^f^^^ + ^^'^^^'^ + 4^^^^^*^- (77) 

The equation for h^^, on brane 1 in the zero mode approximation looks like 



= -|£^(V-^W^)- (7^ 
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Let us consider a static point mass with energy-momentum tensor too = M6{x),tok = 
tij = 0, M denoting the inertial mass again. Then we can find a solution of equation ( |T5D 
for /loo-component 



where r = Thus, we get 



Now let us show that the Randall- Sundrum model in the zero mode approximation is 
equivalent to the linearized Brans-Dicke theory (apparently, it was first noted in Ap- 
propriate formulas for the light defiection in the Brans-Dicke theory are well known and can 
be found, for example, in ||T7|. 



The equation for the fiuctuations of metric in the linearized Brans-Dicke theory looks 

like 

° (V - lv,uh>j = -16nG - (v,. - ^ 

where u is the BD-parameter. It is easy to see that hf^i, satisfies the de Donder gauge 
condition d'^ {h^y — \ri^y}i) = 0. The light defiection angle is given by ]17 



A, . (^) . (82) 



ro \2uo + 4: 

with To being the impact parameter and M being the mass of the static point-like source. 

Now let us examine the four- dimensional effective theory on brane 1. It is easy to see 
that h^^^ (^) does not satisfy the de Donder gauge condition. We can impose it with the 
help of the gauge functions (^). 

In this gauge equation ( |75D takes the form 

□ - lv^..hj = -levrGi (^V " ^ (^m- " ^) ^) • (§3) 

Here and in what follows we drop the superscript of the zero mode approximation. By 
comparing (|83| ) with (|8TD , we get <^ = | (e^'^'^ — l). Substituting this into (|8^), we find the 
angle of light defiection by a static point-like source on brane 1 

^ -r-^-z^ ■ (84) 



The second term in the denominator of formula ( ^4|) and the second term in formula (|80|) 
correspond to the contribution of the radion. One can see that the contribution of the radion 
is e^''^ times smaller than the contribution of the massless graviton. 

Direct calculations lead us to the following results for the remaining cases. All the 
reasonings concerning Galilean coordinates are the same as in section 2. 
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1) We live on brane 1, the mass is located on brane 2 ("shadow" matter) 



□ ( V - 2^A«-^ ) = -levrGie-'^^ (V - 3 - ^ ) M • (85) 



A,.^^^^^4MG.e--f^U. (86) 



The light deflection angle is 

Newton's Law looks like 

^._W-^£^-(l-Hi)G,e-M. (87) 
3 r V 3/ r 

The second term in the denominator of formula (^) and the second term in the brackets in 
formula (|87|) correspond to the contribution of the radion. One can see that the contribution 
of the radion is of the same order, as the contribution of the massless graviton, but the effect 
of the "shadow" matter is e^^ times smaller, than the effect of the ordinary matter (compare 
with (^) and (pUD). It means that the effective theory on brane 1 is phenomenologically 
acceptable, but not interesting (see p). 

2) We live on brane 2, and the mass is located on the same brane 

with Go = Gk- 



g2kR_l • 

The light deflection angle is 



AMG2 , ^ , 



2kR\ 

V = -GAl+^-—]—. (90) 



Newton's Law looks like 



The second term in the denominator of formula ( pUj ) and the second term in the brackets 
in formula (0) correspond to the contribution of the radion. Now the contribution of the 
radion is e^^^ times stronger than the contribution of the massless graviton. It means that, 
in the case of the massless radion, scalar gravity is realized on brane 2. However, if one 
assumes some mechanism for generating the radion mass, for example, the Goldberger-Wise 
mechanism |]19[|, gravity in the zero mode approximation becomes tensor and the hierarchy 



problem solves. 

3) We live on brane 2, the mass is located on brane 1 ("shadow" matter also) 



□ ( V - 2^M-^ ) = -levrGse'^" ( V " 3 [v,. -^]t]. (91) 
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The light deflection angle is 



Newton's Law looks like 



As before, the second term in the denominator of formula ( p2D and the second term in the 
brackets in formula (|93D correspond to the contribution of the radion. We can see that the 
contribution of the radion is of the same order, as contribution of the massless graviton, 
but the contribution of the "shadow" matter to the Newtonian limit is e^^ times smaller, 
than the contribution of the ordinary matter, because of the interaction with the massless 
radion (compare with (|90|)). But if the radion is massive (otherwise we have scalar gravity 
in "our" world on brane 2), the contribution of the "shadow" matter to the Newtonian limit 
is e^^ times stronger, than the contribution of the ordinary matter. It could lead to some 
observable effects, resulting from the distribution of matter in the "mirror" world (brane 1). 

One can shaw that equations (PP]), (P7|), ( PP] ) and ( P5D coincide with the zero-mode parts 
of equations (^2]), (|65D, ( |64D and (|62D respectively. 

We would like to note that in paper [0 the effects of the "shadow" matter were also 
considered. In this paper the deflection angles with the same impact parameters and by the 
same Newtonian masses were compared (Newtonian mass is the coefficient in front of the - 
term in Newton's Law). We have compared the effects of the "shadow" and the ordinary 
matter, as it was made in 0, and we have got the following results: 44% instead of 25% in 
[0 for brane 1 and the difference of the order of e"^^^ instead of 25% in |^ for brane 2. This 
discrepancy appears because we use the Galilean coordinates on brane 2. In the case of the 
massive radion there is no difference in the zero mode approximation at all for both branes. 
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